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Towards a new tool for the indirect detection of Dark Matter : 
building of a SuSy spectrum generator based on micrOMEGAs 

! Pierre Brun^ 

l/^ ' Abstract 

In the quest for indirect signals from dark matter annihilation, powerful computation codes are required. 
I report here a new code based on micrOMEGAs devoted to the analysis of such signals in term of Supersym- 
^ , metry. It computes gamma rays and positrons fluxes in a general SuSy model, as well as the other charged 

' cosmic rays and neutrinos source terms. This work aims to propose an alternative to the DarkSUSY code 

QQ I by providing inclusive signals from SuSy for dark matter indirect searches. Therefore it can be used for 

fT^ ■ sensitivity studies and data analysis. 

m 
o 
\o 
o ^ 

^ ■ 1 Supersymmetric dark matter 

Q ■ 1.1 Supersymmetry and cosmology : the micrOMEGAs code 

I The Standard Model of particle physics describes two distinct families, bosons and fermions as interaction and 

■ matter particles respectively. The integer and half-integer spins can be merged in a unified model if one assumes 

that a symmetry exists between those two sectors, this is the purpose of Supersymmetry (SuSy). In addition to 
this aesthetic issue, SuSy can solve some of the Standard Model limitations, such as the stability of the Higgs 
Kn] ' boson mass as regard to radiative corrections, the unification of forces or the understanding of the electro- weak 

. symmetry breaking mechanism. Since no SuSy particle has ever been detected, it must be that the symmetry is 

broken at today's reached scales, and partners are heavier than the standard particles. Furthermore, assuming 
the conservation of a new quantum number called R-parity, SuSy predicts the existence of relic particles from 
the Big-Bang, these are stable, neutral, weakly interacting, massive and non baryonic (WIMP's). When it has 
these characteristics, the Lightest SuSy Particle (LSP) is a good candidate for dark matter. This possibility 
is very interesting since WMAP measurement of the cosmic microwave background anisotropics point out 
that there must be a large amount of non baryonic dark matter in our Universe. These measurements, in 
addition to the observation of type la supernovae, large scales structures features and constraints from Big- 
Bang nucleosynthesis give the following picture of our Universe content |2]: 

• The universe is globaly flat, with fltot = 1-02 ± 0.02. 

• The major part of the energy is under the form of a cosmological constant n\ = 0.73 ± 0.04. 
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• The matter component is ^matter — 0.27 ± 0.04 and is 84% non baryonic {ilbaryons — 0.044 ± 0.004). 

micrOMEGAs purpose is to accurately compute the rehc abundance in a general SuSy model (31 (namely the 
Minimal Supersymmetric Standard Model, MSSM). To do so, one has to solve the Botzmann equation: 



—— — ~SHn- 
dt 



<av> {n' - 



(1) 



Where n is the number density of relic particles, < av > their annihilation cross section and H the Hubble 
constant. Basically this equation indicates that the annihilation rate has to be balanced against the expansion 
of the Universe, leading to a freeze-out of the LSP. Figure ^ shows how this freeze-out occurs, we can see the 
evolution of the co-moving wimp density. In the first part, the particles interact sufficiently one with another to 
make the LSP annihilate, and the LSP density goes down. At a particular temperature, the Universe expands 
so that the Hubble radius get greater than the mean interaction length. From this moment, the annihilation 
rate is very low and the relic particles density remains constant. 




K=m/T (time 

Figure 1 : Evolution of the co-moving wimp density with temperature 

It happens that some particles have a mass close to the LSP's mass and then contribute to the decrease of 
the LSP density, these processes are the so-called coannihilations 0. In that case one has to replace < av > 
in eq. ^by a sum over all coannihilation channels: 



dn 



N 



(2) 



One great performance of micrOMEGAs is to dynamically include -when they are relevant- all possible coanni- 
hilation channels. The total number of processes that can be involved is about 3000. 



1.2 Supersymmetry breaking 

Once "supersymmetrized" , the Standard Model becomes the Minimal Supersymmetric Standard Model (MSSM). 
if SuSy was effective, Sparticles would have the same mass as the corresponding particles. Superpartners have 
never been observed and it must be that SuSy is a broken symmetry. Then one has to parameterize the break- 
ing by introducing terms in the Lagrangian that explicitly break SuSy. This breaking imposes the introduction 
of over 100 free parameters. In order to make physical predictions, one can fix them or assume a model for 
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SuSy breaking that generates these soft SuSy breaking terms. In this note, we wiU work in two SuSy breaking 
models ft'amcworks : 

• Minimal SUperGRAvity, in which SuSy is broken via gravitationnal interaction. In the mSUGRA model, 
all physical quantities can be derived ft'om 5 parameters specified at the Planck scale, these are : 

— toq: common scalar mass 

— 7711/2 : common fermion mass 

— Aq : universal trilinear couplings 

— tg{f3): ratio of the neutral Higgs vacuum expected values 

— sgn{ii): sign of the Higgs mass parameter fi 

• Amonaly Mediated SuSy Breaking, in which SuSy is broken via the super- Weyl anomaly effects [S] In 
AMSB, 4 Planck scale parameters are enough to describe the MSSM : 

— niQ: common scalar mass 

— gravitino mass 

— tg{P): ratio of the neutral Higgs vacuum expected values 

— sgn{ii): sign of the Higgs mass parameter /i 

In those frameworks, the LSP is the lighest neutralino Xi (or x)i ^ mixing of gauge bosons partners. In 
this note only mSUGRA and AMSB are considered, but the code also allows to work in the Gauge-Mediated 
SuSy Breaking scenario. 

1.3 Cosmological constraints on mSUGRA 

The initial goal of micrOMEGAs is to compute accurately the relic density, and compare it to cosmological 
measurements thus putting constraints on the parameter space of a given model. To illustrate this, a scan 
of the parameter space has been performed. On Figure 12 only the WMAP Aa allowed points are plotted in 
the mo-mi/2 plane, i.e. when the computed thermal rehc density satisfies 0.076 < ilh"^ < 0.148. The mass 
spectrum generator used here is SuSpect |12| . 
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Figure 2: Cosmological constraints on gravity broken SuSy 
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2 AMS02 potential for indirect dark matter searches 



AMS02 is a particle physics spectrometer to be placed on the International Space Station for 3 years. The 
detector allows to measure charged cosmic rays as well as 7 ray fluxes in the range of 1 GeV to a few TeV, 
including particle identification, charge reconstruction, isotopes separation in case of light nuclei W^. It consists 
in different specific sub-detectors, with some redundancy in measurements. A Silicon Tracker surrounded by 
a superconducting magnet provides charge and rigidity, two planes of a time-of-flight detector are used for 
trigger and direction determination. The Transition Radiation Detector and the Electromagnetic Calorimeter 
(Ecal) perform e/p separation and energy determination, and a Cerenkov counter measures charge and /?. 
The principle of indirect search for dark matter is to look for non-standard signals in cosmic spectra, arising 
from annihilations in the local halo 1 . Such a deviation could be seen in the antimatter to matter ratio for 
charged cosmic rays (positrons, antiprotons or antideuteron) , for which AMS capabilities are very high [H]. 
Positron signal is of special interest since HEAT experiment has measured an excess which could be explained 
in terms of a Weakly Interactive Massive Particle (WIMP) annihilation signal 7^. With a background rejection 
of order 10^, an acceptance of 0.04 m? .sr and an energy resolution of 3% in this range, AMS02 is awaited to 
confirm and precisely measure the positron excess In the case of 7 rays, dark matter is expected to produce 
an enhancement of diffuse emission from the halo and possibly point-like sources where its density is high. 
The Galactic center could be such a source, in particular if the halo profile near the supermassive black hole 
is cusped. Another possibility is to observe a GeV scale line emission, which would give a compelling evidence 
of the presence of dark matter since no known astrophysical object could be able to produce it. Although 
originally designed for charged particles detection, AMS02 has high performance in 7 rays detection. Two 
modes allow the spectrometer to perform a measurement whether the photon converts into a e+/e^ pair in 
the upper detector or not. The Tracker is used in the first case and provide good angular resolution (0.05°). 
A specific trigger is used in the other case 9^, for which the Ecal allows a 3% energy resolution and a 1° 
angular resolution. These two complementary ways to detect high energy photons puts AMS02 in great place 
to perform 7 ray astronomy, with high acceptance (of order 0.09 rn? .sr) and a large sky coverage |TP] . 

3 Making micrOMEGAs an event generator 
3.1 Architecture of the code 

The final goal of our package is to predict all indirect signals from SuSy dark matter JI] . The user may choose 
to work in the general MSSM or with any model constrained at the high scale (such as mSUGRA and AMSB), 
in that case the first step is to compute the evolution of physical parameters such as masses and couplings from 
the unification scale to the electroweak scale (EWS) at which cross sections in the halo have to be determined. 
The evolution from the high scale to EWS is managed here by SuSpect |12j which solves the renormalization 
group equations ^. All EWS parameters are given to micrOMEGAs which computes all cross sections. The 
former are discussed in next sections. At this point one has to link these final states to observable particles, 
this task is devoted to PYTHIA |T3j. It computes the partons hadronization- fragmentation and the decays of 
unstable particles. For the 7 signal the integral of the density over a field of view is computed. The charged 
particles have to be propagated through the Galaxy from the location of their production to the Earth vicinity, 
then specific propagation codes are used. In the present note, p and D propagation are not presented, but they 
will be included in the final version. 

^The results presented here make use of SuSpect, but any other RGE code can be used 
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Figure 3: Principle of the code 
3.2 Fluxes from neutralino annihilations 

Dark matter signal in cosmic radiations can be seen either in the matter to antimatter ratio and in 7 rays. In 
this note some results for 7 rays and e+ are presented. Concerning the 7 rays from Galactic center, the flux is 
given by the following expression: 
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l.o.s. 



Two contributions are present in this equation, one purely astrophysical and the other arising from particle 
physics. The right part is the density integral along a line of sight, this depends on the modelling of the dark 
halo. In the other part, from particle physics, the fraction < UannV > I'^n^ contains the total annihilation 
cross section and the neutralino mass, both provided by micrOMEGAs. Here the annihilation cross section has 
to be computed at the temperature of the halo, whereas for the relic density estimation it was computed at 
T freeze-out- dN^ / dE^ is the number of 7 per unit energy for one annihilation and is computed with PYTHIA 
(here version 6.123). Each final state giving different 7 spectra, this is in fact a sum over all final states : 
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(4) 



As X is electrically neutral, it does not couple to photons, therefore there is no tree level direct 7 production. 
Nevertheless, 7 rays arise from the hadronization of the final states partons. A typical example is the decay 
of TT^'s appearing in the hadronization of a quark/antiquark pair. These processes lead to a continuous 7 
spectrum. At the 1 loop level, 7 lines can be produced via processes such as XX ~^ 11 il^^- The 77 and 7Z" 
processes are implemented in our code |15| . and give photons at the following discrete energies: 



77 : i;^, = {N^ = 2) ; 7^° : E^ 



1 



1) 



(5) 



These processes are of special interest since no known astrophysical source could lead to GeV to TeV 7 
lines. The observation of such lines would provide a smoking gun signature for the observation of dark matter 
annihilations in the Milky Way. In this note I present an application of the upcoming tool SloopS developed 
by Boudjema et al. that will compute all of the MSSM processes at the 1 loop level. This package will be 
provided with micrOMEGAs. 
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• Positrons 



For the positron signal, tlie flux is a solution of the propagation equation 



d dnf.+ 
with the source term : 
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(6) 



(7) 



Here, micrOMEGAs provides cross-sections and masses, PYTHIA gives the differential e+ spectrum and a new 
code is used to solve the propagation equation, developed by Salati etal. ^1]. The different parts from particle 
physics and astrophysics are detailed in the following sections. 



3.2.1 Cross sections and final states 

There are 26 possible 2-body tree- level final states for neutralinos self annihilation. These are fermion/antifcrmion 
pairs or allowed combinations of gauge and/or Higgs bosons. Each final state will give its own particles yields, 
with its specific spectral features as we will see in section 4. Therefore it is important to know for a given 
model or set of parameters which final states will be favored. In the mSUGRA framework. Figure 0] shows 
the occurrence of some final states (bb, ti, t^t~ and W^W~) in the too-™i/2 plane with all other parameters 
fixed. 



bb final state probability - A„=0, tg(P)=10, n>0 "c+r final state probability - A„=0, tg(|3)=10, n>0 




m,,, (GeV) m„, (GcV) 



Figure 4: Final state occurrences with Aq = 0, /i > and tg{(3) — 10 in mSUGRA 
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It appears from Figure 0] that a typical mSUGRA case is the one in which annihilation into a bb pair is 
dominant, followed by the annihilation into r's. Some regions of the parameter space lead to an annihilation 
into top quarks, with a threshold effect, or into gauge bosons such as W^s. However, by comparing this figure to 
Figure 12 one can see that the annihilation into W^W" is not favored by cosmology in the case of mSUGRA. 
However it will be shown in last section that other supersymmetry breaking scenarios can lead to this final 
state. 

3.2.2 Dark halo modelling 

The exact local dark matter density is not known and one has to parameterize the halo density profile in order to 
match the observed gravitational effects. The local dark matter density is in the range pq = 0.2 — 0.8 GeV.cm~^ 
and that does not give information about the density far from the center or near the center. Therefore one has 
to assume a halo profile that gives the proper local density. The following parametrization is implemented in 
the code, it allows to describe number of halo profiles (motivated by numerical simulations). 



PCDKlir) = P0 
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In addition to this, any kind of central shape can be specified by hand, this is of special importance if one 
wishes to test different accretion models onto the central supermassive black hole. Figures [S] and El give the 
most popular halo profiles parameters and shapes, in the case of Navarro- Frenk- White (NFW), Moore and 
isothermal sphere (see and references therein). 
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Figure 5: Halo parameters for three profile types 
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Figure 6: Three examples of halo profiles 



This parametrization allows to predict diffuse fluxes from the halo as well as fluxes from the Galactic center. 
In order to make a prediction for a 7 flux from the Galactic center, one has to integrate the signal inside some 
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solid angle around the direction of the Galactic center. To do so, the code sums all lines of sight contributions 
such as the ones shown in Figured 




Figure 7: Dark matter density along lines of sight with different angles toward the Galactic center 

As an example, the table of Figure El show the normalized value of the astrophysical factor < J >, defined 
as the mean value of J over a given Afi solid angle. Here < J > is : 

< J >An= J Pl{r{h))dl^ (9) 
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Figure 8: Astrophysical terms calculation with the package 



These figures are obtained with our package, they are in perfect agreement with previous calculations 0. 
These figures show that in the case of this observation, a better angular resolution lead to a better observation 
only in the case of a cuspy halo profile. Indeed for a halo with a flat dark matter distribution in its center, the 
mean value of the astrophysical factor does not vary when the solid angle is changed. 



4 First example : 7 rays from Galacic center in mSUGRA 

I wish here to illustrate the use of our package for computing 7 spectra from the Galactic center. The 
astrophysical part used here is a NFW parametrization (the distance to the Galactic center is tq = 8 kpc and 
the local dark matter density is taken at pq = 0.3 GeV.cm~^) . The annihilation rates are integrated inside a 
1° opening angle cone around the Galactic Center. Then the 7 ray flux is given by the following formula : 
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As examples, I use two sets of mSUGRA with parameters detailed below. Both points match relic density and 
accelerator constraints. The first one is used by W. de Boer et al. |THJ and has the advantage of giving quite 
a high signal, with a total cross section of order 10~^^ cm^s~^ . The second point gives less flux but a visible 
line (C benchmark point ^2]). Figure El gives some characteristics of the chosen sets of parameters. 



Set 1 


Set 2 


mo = rai/2 = Aq = 500 GeV 
tg{l3) 50, ^ > 


mo = 113 GeV, mi/2 = 375 GeV, Aq ^ 
tg{P) = 20, fi>0 


n^K^ = 0.098 
mx = 206.9 GeV 

< atotv >= 1.81 .10-26 cTO^s-i 

< a{hh)v 0.87 x < atotv > 

< cr{T+T^)v >= 0.13 X < atotV > 

< a{-/-f)v 2.10"^ X < atotv > 
< a{-fZ")v >= 4.10-6 X < atotV > 


n^h'^ = 0.128 
m^ = 151.5 GeF 
< atotv >= 6.96 .10-28 cm^s-i 

< a{bb) >= 0.58 x < atotV > 

< o-(r+r-) >= 0.39 x < atotV > 

< <^ill) >= 0.01 X < atotv > 

< cT{-fZ") >= 0.001 X < atotv > 



Figure 9: Main features of the two sets of mSUGRA parameters considered in this part 



All the final states cross sections being determined, the 7 spectra for each one of them has to be computed. 
This is done by using micrOMEGAs final states as an input to PYTHIA, which hadronizes the partons, simulates 
the jets fragmentation and the decay of all unstable particles that may appear in the hadronization process. 
Different spectra are obtained for each final states. The final spectra results in weighting these ones with the 
computed cross sections. The user has different ways to obtain the final signal: 

• Monte Carlo generator: The spectrum is produced event by event with a final state random choice 
according to their probability. 

• Weighted channels: Channels are produced separately and then weighted as regard to their cross section. 

• Use of tables : A catalog of spectra is available, those being interpolated to match the specified parameters. 
Although not very accurate, this method is much faster and allows a fast rough determination of the 
spectra, which is useful for scans of the parameter space. 

In this paper the second method is used and all channels with probability greater than 0.1% plus 77 and 7^" 
have been considered in the signal prediction. 
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Set 1 



For the first set of parameters, the following spectrum for the 7 rays arising from Galactic center is obtained. 
Here the flux is high and the 7 lines are undistinguishable from the continuous hadronization spectrum. 




Figure 10: 7 flux from Galactic center for set 1 



Set 2 

For the second set of parameters, Figure [TTI is obtained. As expected the flux is quite low but one can see very 
clearly the lines induced by xx ~* 77 ^-nd XX ~* 7^*^- 
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Figure 11: 7 flux from Galactic center for set 2 

Of course, those two spectra have to be added to the background and compared to AMS02 acceptance. 
According to the background expectation, it will be less obvious to extract the signal in the latter case (in next 
section a estimation of AMS02 sensitivity to this signal is presented) . One can notice that the specific spectral 
features can only be observed with a very good energy resolution. In order to illustrate this, one can plot 
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the same spectrum with an typical Atmospheric Cerenkov Telescope (ACT) energy resolution, as it is done in 
Figure ^1 (a ~30% energy resolution is assumed) . 



In the case of an ACT, the flux suppression at could be seen, but the 7 line itself is too much smeared 
by the energy resolution. Only a space mission like AMS02 would be able to do so. 

Estimation of theses fluxes observability with AMS02 

As it was said before, the first set of parameters provides quite a high signal. In order to have a rough estimate 
of the observability of this flux with AMS02, an average acceptance of 0.09 m^.sr is assumed, and the signal 
is integrated inside a solid angle of 10~^ sr corresponding to the 1° angular resolution of the Ecal. For 3 years 
of data taking, the measurement of the 7 flux in case of set 1 would lead to the points of Figure E| Here the 
expected background is a power law of spectral index 2.72 extrapolated out of the EGRET measurement of 7 
rays from the Galactic center region P^ . 
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Figure 12: Same as previous Figure ITTl with an ACT energy resolution 
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Figure 13: AMS02 measurement of the flux from Galactic center for set 1 
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In the second scenario, the signal shown on Figure ^2 gives a signal that leads to an excess at the level of 
less than la over the background, even in 10 years of data taking. It would be possible to observe this line if 
the dark matter halo were more cusped than a typical NFW profile. In Figure the observation of the line 
by AMS02 in 3 years is presented under the more optimistic assumption of a central cusp of index 7 ~ 1.25. 
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Figure 14: AMS02 measurement of the flux from Galactic center for set 2, with a speciflc halo proflle 

5 Second example : positron signal in AMSB 

In the AMSB scenario, the neutralino is a pure wino, and the annihilation cross section is always higher than 
in the mSUGRA case. This has two main implications : 

• The thermal relic density is very small so that one has to assume non thermal x production. Such a 
non-thermal relic density can be caused by non-standard cosmology |19| or decays of gravitinos produced 
at the end of inflation jSHj. 

• The XX ~^ W^W^ channel is always dominant as soon as it is kinematically allowed. 

The annihilation into W bosons gives harder e"*" signal than in the typical bb case of mSUGRA. That is 
why this framework has been chosen to illustrate the generation of positron signal. In the example here, the 
SPS9 benchmark point (21jis used, for which : 



SPS 9 



mp ^ 450 GeV, TO3/2 = 60 TeV, tg{P) = 10, /i > 



(fl^h )thermal — 0.0018 

= 197.6 GeV 
< cr^^tw 1.93 .10-2"* cm^s-i 
< a{W+W-)v >= 0.991 X < cjtotv > 
< a{jj)v >= 0.0021 X < atotv > 
< (T{"fZ°)v >= 0.0065 X < atotv > 



Figure 15: Main features of the SPS9 set of parameters 



For the flux calculation, the same halo model as above is used. Figure ^| shows the e"*" before propagation 
and after propagation. The code that was used here is currently under development 14^ and its final version 
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will be implemented in micrOMEGAs. It makes use of a new positron propagator for the diffusion in the galaxy 
and its spatial resolution should not be limited, in order to allow the study of small dark matter structures. 
However for the signal shown here, the e+ source term is integrated over the whole galactic halo, taken to be 
smooth. 




Figure 16: Positron source term and flux after propagation (see text) 



Figure 1161 shows the superimposition of the e"^ source term (blue curve) and the flux after propagation. 
Here the red curve has been re-scaled in order to see the spectrum distortion due to propagation. The scale on 
the left corresponds to the source term in units positrons.GeV~^ . The scaling factor used here to draw the flux 
curve is of order 10^ crn?.s.sr. Dark matter substructures (clumps) are awaited to produce an enhancement 
of the signal ^1]. Because of the high annihilation cross-section and the specific high energy features of the 
W-induced positrons, a boost factor of only order unity is required to fit the HEAT positrons excess [3 . Figure 
II 71 shows the SuSy signal in term of positron fraction, compared to the HEAT excess. The expected background 
is taken from |22|. 




Figure 17: Fit to the HEAT data in AMSB 
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In Figure nn the supersymmetric signal is multiplied by a factor of 3.8, this figure is to be compared to the 
typical mSUGRA case, in which the boost factor can vary from 10 to over 1000. An important feature of the 
AMSB scenario is the charginos/neutralinos mass degeneracy, which implies a very low sensitivity for hadron 
colliders to this model. Therefore if AMSB is the proper description of high energy physics, AMS02 would be 
in a very good position to observe SuSy. 



6 Other charged channels : antiprotons and antideuterons 

For the moment, the propagation of p and D is not implemented in the code. However, the injection spectra 
that stands for the source terms in the propagation equations are available. Figure IT^ shows the source terms 
for these two channels. 




Figure 18: p and D source spectra 

The D spectrum is computed according to the coalescence model, see ref [^S]. The adjunction of these 
channels in the code is currently under way. 

7 Conclusions and outlook 

In this note I gave some examples of SuSy signals which could be seen by AMS02. These illustrations show the 
potential of the package we are developing. Its possibilities will be even greater than what is shown here, it 
allows to work in any SuSy model and with any halo profile. In the final version, substructures and extragalactic 
sources studies will also be possible. Some comparisons with the only existing similar tool DarkSUSY |E3 are 
being performed and we are currently working on the implementation of the p and D propagation. In the 
near future, the nature of dark matter could be unveiled by the spectral features of the annihilation signals 
in the Milky Way. Once published, this package could be a powerful tool for data analysis of all experiments 
performing indirect searches of dark matter. 
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